Weak Lensing Science, Surveys, and Systematics 
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Weak gravitational lensing is one of the key probes of the cosmological model, dark energy, and 
dark matter, providing insight into both the cosmic expansion history and large scale structure 
growth history. Taking into account a broad spectrum of physics affecting growth - dynamical 
dark energy, extended gravity, neutrino masses, and spatial curvature - we analyze the cosmolog- 
ical constraints. Similarly we consider the effects of a range of systematic uncertainties, in shear 
measurement, photometric redshifts, and the nonlinear power spectrum, on cosmological parameter 
extraction. We also investigate, and provide fitting formulas for, the influence of survey parameters 
such as redshift depth, galaxy number densities, and sky area. Finally, we examine the robustness 
of results for different fiducial cosmologies. 
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I. INTRODUCTION 

Weak gravitational lensing is a key probe of cosmology, 
in terms of both the properties of spacetime as a whole 
and the distribution of matter within the universe P, ■ 
It has the characteristic of being sensitive to both the 
cosmic expansion and the growth of structure, making 
it ideal for testing the consistency of these two windows 
on the universe, e.g. as predicted by general relativity. 
Moreover, weak lensing has substantial complementar- 
ity with other probes, enabling more precise constraints 
through breaking degeneracies. 

Because of this richness, weak lensing not only probes 
many aspects of physics but is also affected by them. Ne- 
glecting some of the sources impacting growth of struc- 
ture could bias the results, or at best underestimate the 
true uncertainties. In this article we present a series of 
analyses calculating the effects on weak lensing science 
when incorporating simultaneously the major possibili- 
ties for modification of ACDM growth. These include 
using dynamical dark energy, extended gravity, neutrino 
masses, and spatial curvature. This generalizes investi- 
gations that only consider one or two of these physical 
ingredients at a time. 

We also examine how conclusions on survey properties 
alter when the fiducial model is different from ACDM, 
and in the presence of systematics in the photometric red- 
shifts, shear measurement, and nonlinear power spectrum 
form. In Sec. |ll] we present our methodology, includ- 
ing the basic equation for the weak lensing shear power 
spectrum and some hidden assumptions in its derivation, 
the set of cosmological parameters, and the framework 
of our Fisher analysis code. Section IIIII outlines the 
survey parameters and systematic models, including a 
parametrized form for the nonlinear matter power spec- 
trum. Cosmological parameter constraints are discussed 
in Sec. lIVl including effects of an extended parameter set, 
the fiducial cosmology, survey parameters, and systemat- 



ics. We also present a selection of trade studies and the 
resulting fitting functions, that could guide more effective 
survey design. 



II. FOUNDATIONS AND METHODOLOGY 
A. Weak Lensing Shear 

We briefiy review the key equations for the weak lens- 
ing shear power spectrum. For a comprehensive treat- 
ment, see The shear (equivalently the convergence 
in the weak lensing limit) is a weighted integral of the 
mass density field. The average convergence of a light 
ray bundle from sources in a source bin i, in the sky 
direction 6, is 



dr]W.,{ri)S{ri,rie) , 



(1) 



where S ~ 5p/ p is the fractional mass density perturba- 
tion and rj is the conformal time (the upper limit will be 
cut off by the kernel). The convergence weight function 
or kernel is 



W,{z) = -nr„Hl{l + z)r^{z) 



'rf.'^^. (2) 



Here the redshift z is another way of measuring the con- 
formal distance rj, r]{z') is the distance to the source, 
7y(z, z') the distance between lens and source and ri{z) the 
distance to the lens. The number of sheared objects (in 
source bin i) along the line of sight per unit sky area and 
per unit redshift is given by the source distribution ni {z), 
with the total number per unit sky area nf = J dz ni{z). 
Finally, the matter density in units of the critical density 
is f2,„, and Hq is the Hubble constant. 

By Fourier transforming to multipole space and form- 
ing the two point correlation, one obtains the power spec- 



2 



trum (omitting source bin subscripts for now): 

= i2nfS^ir+e')Ce, (4) 
in analogy to the mass power spectrum Pk, 

{5{k)5{k')) = {2nf5°{k + k')Pk. (5) 

Note that 5^ denotes a Dirac delta function and not a 
density perturbation and we employ the flat sky approx- 
imation. 

From the above equations one can see that, in the Lim- 
ber approximation, the shear power spectrum between 
two samples (source redshift bins) i and j is 

Cl^j^ I ^W,{z)Wj{z)Pk{k^i/7j;z). (6) 

One immediately sees that weak lensing is sensitive to 
several distinct quantities: cosmic distances in terms of 
r]{z) and r](z, z') and hence the expansion, the present 
physical mass density D.mH'^, the matter power spectrum 
Pk{z), which involves the expansion history $7^(0), the 
growth history, and the linear to nonlinear density map- 
ping, and the source distribution ni{z), which will re- 
quire translating photometric redshifts into true ("spec- 
troscopic" redshifts), and allows the use of redshift to- 
mography or crosscorrclations between distinct samples 
at different stages of growth. 

A number of assumptions are implicit, and often for- 
gotten, in the derivation of Eq. ([5]). To obtain the rela- 
tion between the convergence and the matter density per- 
turbation in Eq. (P) one solves the geodesic equation in 
general relativity specialized to a Friedmann-Robertson- 
Walker (FRW) spacetime. A change in the gravitational 
light deflection law, such as through unequal time-time 
and space-space metric potentials, will alter the result, 
as will a modification of the Poisson equation relating 
the potential to the density perturbation. Properly, the 
kernel W should be multiplied by the G function gener- 
alizing Newton's constant (sec j5[). See @ for more de- 
tails; here we assume that both relations take their stan- 
dard form. Likewise we assume that only matter clusters 
significantly. One could also explore the effects of tak- 
ing the light geodesic in a slightly inhomogeneous rather 
than FRW spacetime, i.e. one where the light rays re- 
ceived by the observer preferentially traverse underdense 
regions. This then changes the distance factors from the 
FRW to the Dyer-Rocder form 0. While this is for- 
mally required for consistency, practically it is insignifi- 
cant for weak lensing @ . Finally, in the conversion from 
the homogeneous mass density to flmHQ one assumes 
Pmia) ~ o,~^j that is conservation of matter (no interac- 
tions with dark energy, for example). To summarize, the 
standard formula must be corrected when investigating 
cosmological models involving certain types of gravity, an 
inhomogeneous universe, or matter interactions. We will 
keep all the standard assumptions. 



B. Physical Ingredients 

In weak lensing the cosmological model enters through 
both the expansion history (distances) and the growth 
history (matter power spectrum Pk). The cosmological 
parameter estimation must therefore take into account all 
parameters that could impact these. On the expansion 
side this includes not only the matter density but dy- 
namical dark energy through its time varying equation 
of state, and the possibility of spatial curvature. On the 
growth side in addition to the usual scalar perturbation 
power index one must consider again the effects of dy- 
namical dark energy and spatial curvature, plus neutrino 
mass and the possibility of growth modifications through 
extensions to general relativity. 

Thus the vanilla ACDM cosmological parameter set 
must be enlarged to account for these physical effects. 
We include the dynamical dark energy equation of state 
w{a) through the two parameters wq and Wa, where 
w{a) ~ ti;o+Wa(l— a) where a = 1/ (I +z), the spatial cur- 
vature density flk , the sum of neutrino masses ^ , and 
modified growth through the gravitational growth index 
7, especially suitable for gravitational modifications that 
are scale independent on the scales relevant for weak lens- 
ing. Dark energy suppresses growth due to the increased 
Hubble expansion rate and smooth spatial distribution 
(we do include dark energy perturbations but this con- 
tributes little to the matter power spectrum). Neutrino 
mass suppresses growth through free streaming. Spatial 
curvature acts as an unclustered component and so effec- 
tively dilutes the matter clustering. Gravitational modi- 
fications can enhance or suppress growth. Cosmological 
parameter estimation when including only one of these 
effects can lead to incorrect conclusions, if others exist as 
well, due to the similarity of their physical impacts. 



C. Weak Lensing Code 

The weak lensing cosmology code developed by the au- 
thors is possibly unique in that it includes all four of the 
cosmological infiucnccs on growth, together with inclu- 
sion of observational systematics. This paper presents 
examples of analysis using this Berkeley weak lensing 
code to study the impact of these additional contribu- 
tions. The code uses a modified version of CAMB 
to include the additional parameters in calculating the 
linear matter power spectrum, Halofit generalized by the 
inclusion of the growth index 7 to form the nonlinear 
matter power spectrum, and then integrates over the ge- 
ometric and source distribution kernel to compute the 
weak lensing shear power spectrum. Cosmological pa- 
rameter estimation is carried out through Fisher matrix 
analysis, with the error covariance matrix including sys- 
tematics in photometric redshifts of the source galaxies 
being Icnscd and in the shear measurements. We also 
consider systematics in the form of the nonlinear mass 
power spectrum in Sec. IIII CI 
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TABLE I: Parameter Set and Fiducial Values 



Cosmological 


O y oLcllicLLiL- 


Parameter Value 


Parameter Value 




Photometric 


iZj/l U.U/200 
fic/i^ 0.1109 

fide 0.734 
Us 0.96 
as 0.8 


{0.0, 0.3, 0.6, 0.9, 1.2, 1.5, 1.8, 2.1, 2.4, 2.7, 3.0}'^ 
scatter, crj 0.03 (1 + ^z) 
bias, Zbias 0.0 (1 ZpJ 


Additive Shear' (5 tomographic bins) 


a' 1.0 
p 0.5 
^. 1000 

b' {lo~^ lo-^ lo-^ lo~^ lo-""'} 


Beyond-ACDM 


0.15 eV^ 
Qk 0.0 
7 0.55 
Wo -1.0 

0.0"= 


Multiplicative Shear'' (5 tomographic bins) 


fi {0.0, 0.0, 0.0, 0.0, 0.0} 



"The photo-z systematics are defined through linear interpolation between these 11 redshift nodes. 

'Defined through Eq. ll23l l. 

■^Three degenerate species assumed. 

"^Defined through Eq. ([TSJ. 

'^When WQ and Wa are varied, their fiducial values are slightly offset so as to prevent w{a) crossing — 1. 



The Berkeley weak lensing code will be made publicly 
available in the near future. 

The cosmological parameter set for the Fisher matrix 
analysis includes the vanilla ACDM parameters, namely, 
the density parameter for baryons f2h/i^, for dark matter 
flch^, and for dark energy ^ldc^ a-s well as the scalar spec- 
tral index Ug and the amplitude of matter fluctuations at 
redshift z = on 8 h~^Mpc scales, erg- Beyond-ACDM 
parameters explored in this paper include the growth in- 
dex 7, the sum of neutrino masses ^ m,^, spatial curva- 
ture parameter Qk , and the dark energy equation-of-state 
parameters wo and Wa- Although not discussed in this 
paper, the code also includes the following parameters: 
the number of relativistic species A^eff, reionization opti- 
cal depth T, and the running of the spectral index a, in 
anticipation of synergistic studies with CMB and CMB 
lensing data sets. 

Apart from the cosmological parameters, a set of pa- 
rameters describing systematic effects can be varied. 
These are further discussed in Sec. IIIII Table U dis- 
plays the cosmological and systematic parameters and 
their fiducial values. 

Information on parameters {p^} can be obtained, un- 
der the assumptions of Gaussianity of the observables, 
from the Fisher matrix 

where to simplify notation we define 

X£,a=z{i-l)/2+] = Ci.ij {i > j) ■ (8) 



The covariance matrix is given by 

{Cov()ab = (2^+ l)f^ + C^iiC^/fc] ' (9) 

with a = i{i - 1) /2 + j and b = k{k - 1) /2 + I. The total 
spectra C^"' arc defined as the sum of the cosmological 
signal and noise: 

Ce% = Ce,,, + N,,,j (10) 

where the noise spectra Ni encode the effects of shape 
noise and other systematics, as discussed in Sec. IIIII For 
the purpose of this study, we have set £niin = 2 and 
^max = 3000, the latter being the maximum multipole 
up to which the Gaussian approximation is sufficiently 
valid. Here /sky is the fraction of sky covered by the ex- 
periment; we chose /sky = 0.097 corresponding to 4000 
square degrees. 



III. SYSTEMATICS 

The number of galaxies whose shears are measured in 
a large survey can be many millions so the statistics are 
copious and systematic uncertainties play a correspond- 
ingly larger role. We do not present a comprehensive 
analysis (see, for example, [13] ) - many of the effects de- 
pend on observational details - but include the two stan- 
dard sources of uncertainty coming from the data and a 
further one from theory. 
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A. Photometric Redshift Systematics 



B. Shear Systematics 



Obtaining accurate, spectroscopic redshifts for the 
many millions of measured galaxies is impractical, so 
weak lensing surveys rely on photometric redshifts. 
These will be imperfectly calibrated by a spectroscopic 
subsample, and have a residual scatter and bias with re- 
spect to the true redshifts. Since this alters the kernel in 
the shear power spectrum, these systematics will propa- 
gate into the cosmological parameter estimation. 

We follow the photometric redshift model as described 
in For an observed photometric redshift Zph, the 

probability distribution p(zpii\z) of the true redshift is 
modeled as a Gaussian distribution 



PiZphl 



1 



■ exp 



{z - Zph - Zbias)^ 



2cr? 



(11) 



where zua.siz) and crz{z) are the bias and scatter in 
p{zpii\z), respectively, and where we allow them to be 
an arbitrary function of redshift z. 

We then assume an overall true galaxy redshift distri- 
bution n{z) = (PN/dzdVL of 



n{z) (X z"exp [— (z/zq)'^ 



(12) 



where we have implemented a = 2 and (3 = 1; zq is 
a characteristic redshift which is related to the median 
redshift by zq = Zmed/2.674 in that case, and the nor- 
malization is fixed by the total number of galaxies per 
steradian 



n" = / dz n{z). 

'0 



(13) 



The true distribution for objects binned in photometric 

redshifts, with lower to upper limits of z^y^ to z^^ , is 
then 



ni{z) 



,(« + !) 
-ph 



dzphn(z)p(zph|z) 



(14) 



= -n(z)[erf(x,+i) - erf(a;j)] (15) 



where Xi is defined as 



-ph 



Z + Zbi 



s) /V2<j, 



(16) 



The total number of galaxies per bin per steradian is then 
.A_ I (17) 



nj' = / dz ni{z). 
'o 



We identify Zbias and as the nuisance parameters for 
the Fisher matrix analysis, where they are taken to be 
z dependent. These values are defined at npzbin inter- 
vals, where intermediate values are linearly interpolated. 
Hence there are 2 x npzbin nuisance parameters associ- 
ated with photometric redshift-related systematics. It 
is the uncertainties in these parameters, more than the 
values of the parameters per se, that impact the cosmo- 
logical parameter estimation. 



Measurement of the shear from galaxy images is a com- 
plicated process, since galaxies have (unknown) intrinsic 
shape cUipticities and imperfect resolution and optical 
distortions due to the telescope, detectors, and atmo- 
sphere contribute as well. The galaxy shapes can be 
treated statistically through a shape noise contribution 
cr^ / nf in Eq. pUj) but the other uncertainties give resid- 
ual systematics. 

We follow the formalism in [l^ for additive and mul- 
tiplicative shear systematics in the weak lensing signal 
extracted from the galaxy images. The multiplicative 
systematic is typically generated from the incomplete re- 
moval of the finite-size effects of the point spread function 
(PSF), but can also be generated from non- ideal weight- 
ing of galaxy shapes to estimate the net shear. Its effect 
can be encapsulated in the factor / 



7(z,n) = 7(2:, n) [1 + f{z,n) 



(18) 



where 7(2, n) and 7(2, n) are the estimated and true 
shear at some (true) redshift z and direction n. Notice 
that / is itself is a function of z and n; it is also time 
dependent (due to the PSF dependence on these quan- 
tities). Within the tomographic bin i, we average the 
multiplicative factor / over n and z such that we are left 
with the bias fi = (/(z,;,n)). The observed shear corre- 
lation function between the ith and jth tomographic bin 
is then 

(7(zi,n)7(zj,n + dn)) (19) 
= (7(zi, n) 7(zj, n + dn)){l + /,)(! + /,) 
- (7(^i,n)7(2j,n + dn))(l + /, + /j). 

The additive systematic is typically generated from the 
anisotropy of the PSF, and is characterized by 7add, 



7(z, n) = 7(z, n) + ji,dd{z, n) . 



(20) 



Since the PSF anisotropy is uncorrclatcd with the true 
shear 7, we can assume (7(n)7add(n -t- dn)) = 0. The 
non- vanishing term is then (7add(n) 7add(n + dn)), with 
Legendre transform -Padd(^)- The additive shear system- 
atic then affects the weak lensing angular power spectra 
as 



add • 



(21) 



It is the power spectrum of the residual after character- 
ization of the additive shear that impacts the observa- 
tions. Hence the additive shear error is given by 



7add(2i,n) = hir{n) , 



(22) 



where r(n) is a direction dependent random variable, and 
hi is the characteristic additive shear residual amplitude. 
The additive shear error power spectrum is assumed to 
take the form of a power law [l^ , 



(23) 
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where the coefRcient p is set to unity for i = j, and 
fixed to some fiducial value for i ^ j- The choice of 
is arbitrary and is degenerate with the parameters 
bi; we choose = 1000 since the weak lensing signal is 
generally sensitive to this region. We take a' = 1; [l2j 
find their results are insensitive to this power law index. 

We include fi, bi, and a' as nuisance parameters for 
the Fisher matrix analysis. The parameters fi and bi 
are defined for each tomographic bin. Hence there are 
2xntom+l nuisance parameters associated with shear 
systematics, where ntom is the number of tomographic 
bins. 



C. Nonlinear Power Spectrum 

Weak lensing accesses information from both the lin- 
ear and nonlinear matter density regimes. The nonlinear 
matter density power spectrum is a key element in the 
shear power spectrum, but the contribution of the non- 
linearities is imperfectly known. Frequently the nonlin- 
ear power spectrum is generated from the linear power 
spectrum through the Halofit prescription p^ . but even 
for ACDM this is accurate at only the ~ 10% level. A 
detailed and valuable investigation of the effects of sys- 
tematic uncertainties in the nonlinear power spectrum 
on cosmological parameter estimation from weak lensing 
was carried out by [3] . 

They considered the density power spectrum as piece- 
wise constant in bins of wavenumber k and found that 
each bin needed to be constrained to ~ 1% level in order 
not to degrade cosmology estimation from next genera- 
tion weak lensing surveys. The worst bias caused by the 
systematic uncertainty comes from step features in the 
power spectrum, i.e. sharp changes with k. Many physi- 
cal infiuences on the power spectrum, however, would be 
expected to be smooth with wavenumber. 

Here, we consider one specific case of a smooth mul- 
tiplicative function representing the systematic uncer- 
tainty. The form adopted is roughly motivated by mod- 
ified gravity effects on the matter spectrum as found 
from both analytic arguments and numerical simulations 
[ist . Although inspired by f{R) gravity models with 
a chameleon mechanism restoring general relativity on 
small scales, the form should simply be taken for itself, 
representing a deviation from the ACDM prediction of 
Halofit on small scales, then a constant offset on even 
smaller scales out to the fcmax limit included in the weak 
lensing calculations. (There might be a restoration to 
zero offset on even smaller scales past our limit.) 

The example parametrization is 



Pik) = PHF(fc) 1 



1 + {k/kof 



(24) 



where ko = Ih Mpc ^ , Anl is a free parameter and Prf 
is the non linear power spectrum as given by Halofit [l^ . 



The power spectrum systematic enters through assum- 
ing the wrong value of Anl, e.g. using Halofit only, 
that has Anl = 0, when there is really a modifica- 
tion. We first consider the systematic as a bias JAnl = 
^NLassumod " ^NLtrue- This propagatcs iuto Si bias in 
the cosmological parameters through the Fisher bias for- 
mula. In general, if we have A'^-l-1 parameters {pi} and we 
misestimate the parameter p^v+i by dpN+i, the resulting 
bias in pi is 



dpi 



N 



dpN 



+1 



(25) 



In our case 5pN+i ~ <5Anl- 

For our fiducial weak lensing survey, combined with 
Planck CMB information, the multiplicative leverage of 
a power spectrum systematic is a factor of 18 (6) for Wa 
(wo)- That is, a 10% misestimation of ^nl biases Wa by 

1.8(7. 

Alternatively, if one is sure of the form for the 
power spectrum, one can try to fit for Anl- In this case, 
the other parameters will not be biased, but their errors 
might degrade. In fact, we find that for the weak lensing 
plus Planck CMB combination, the degradation is less 
than 5% in the dark energy parameters and the power 
spectrum uncertainty can be determined to (t(Anl) = 
0.016. This highlights the crucial role of understanding 
the form of the k dependence of the nonlinear matter 
power spectrum. 



IV. COSMOLOGY CONSTRAINTS AND 
SURVEY DESIGN 

A. Effect of Systematics 

In this section, we investigate the effect of the various 
systematics discussed in the previous section on the cos- 
mology constraints. We consider a fiducial survey with 
a sky coverage of 4000 deg^ and an effective projected 
galaxy number density rieS = 55arcmin~^. We take the 
intrinsic galaxy shear to be cr^ = 0.27. We assume a red- 
shift distribution given by Eq. p2p with Zmod = 1 and 
subdivide the galaxy sample into five rcdshift bins, given 
by z = 0.6-1.0-1.5-2.1-3.9. It was shown in [ul that 
five tomographic bins is sufficient to capture the rcdshift 
information. 

For the photometric redshift errors (see Sec. IIII A]) , we 
assume a fiducial redshift scatter [ill 



cj^{z) = 0.03(l + z) 



(26) 



and our fiducial bias is zero. To characterize the uncer- 
tainty in the scatter and bias, we parametrize az{z) and 
Zhiasiz) by their values in npzbin= 11 bins, include these 
22 parameters in the Fisher matrix, apply a prior of 0.01 
on each, and finally marginalize over these parameters. 

For the additive shear bias (see Sec. IIII B| . we assume 
fiducial values bi ~ 10^^ for each tomographic bin. Since 
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the true additive shear is unknown (if it were known, 
we would not consider it a systematic), we again add 
the bi as parameters in the Fisher matrix, with priors 
of 10~^. Finally, we consider the multiplicative shear 
systematic (see Sec lIIIBj) . with fiducial parameters fi = 
and priors of 0.01. 

In Fig. [U we show the effects of the systematics on 
constraints for the dark energy parameters wq and Wa, 
marginalized over all other parameters. We also include 
a forecasted Planck CMB prior in the bottom panel only, 
showing how that information breaks certain degenera- 
cies (mostly between parameters other than wq and Wa, 
but this results in narrowing the WQ-Wa contours). 

In either case, the systematics have noticeable effects, 
even for a 4000 square degree weak lensing survey, on 
the dark energy constraints. We show the effects of each 
systematic individually, as well as all of them included si- 
multaneously. The strongest impact comes from the mul- 
tiplicative shear and the photometric redshift bias uncer- 
tainty, at least for the prior levels we have considered as 
reasonable expectations (although one should be able to 
improve on Zbias with a good spectroscopic calibration set 
for the photometric redshifts). In general, we expect the 
systematics contribution is largest from the multiplica- 
tive shear error and redshift bias since these components 
most directly affect the calibration of the lensing signal, 
and hence cosmological parameters. The results agree 
qualitatively with jllL Il2l | , so we now proceed to study 
the extension to a larger array of cosmological physics 
affecting large scale structure growth. 




Weak Lensing only 



— no systematics 




Weak Lensing plus Planck prior 



B. Effect of Cosmological Physics 

As mentioned in the introduction, several distinct 
types of physics affect cosmic growth. In this section 
we investigate the effects of beyond-ACDM parameters 
on dark energy constraints. Besides the dynamical dark 
energy equation of state itself, we consider three param- 
eters describing these influences — the spatial curvature 
density parameter fi^, the sum of the masses of neutri- 
nos ^m,y, and the gravitational growth index 7. The 
effects of varying each of these individually and all at the 
same time are shown in Fig. O Allowing for variation of 
all three parameters together leads to a worsening of the 
dark energy figure of merit (FOM) by a factor of ^ 2. 

Either relaxing flatness (allowing for $7^ 7^ 0) or that 
growth must follow expansion (allowing for 7 7^ 7gr) 
have large impacts. For curvature, recall that this en- 
ters into both the distances and growth part of the weak 
lensing signal; generally combining accurate distance in- 
formation, from a supernova or baryon acoustic oscilla- 
tion survey, would ameliorate its impact on the cosmol- 
ogy constraints. Allowing for growth to deviate from 
the general relativity prediction (with matter being the 
dominant clustering component), weakens the dark en- 
ergy constraints, but guards against substantial bias if 
such deviations do exist [l6| . For example, a departure 



FIG. 1: Dependence of the dark energy equation of state 
constraints (68% CL contours) on the following systematics: 
uncertainty in the photometric redshift scatter and bias, mul- 
tiplicative shear and additive shear. We show contours where 
each effect is taken into account individually, a contour with- 
out systematics, and one with all four systematics included. 
For the priors chosen here, the multiplicative shear and the 
photometric redshift bias are the dominant systematics. Tak- 
ing into account all of them, errors are increased by about a 
factor of 2, and the contour area ( "figure of merit" ) by about 
a factor of 5. In the bottom panel only, we also include our 
code-generated forecasted Planck CMB Fisher matrix. 



A7 can give rise to a bias Awa ~ 8A7. Although this ex- 
tra freedom generically increases the wo-Wa contour area 
by factor of ^ 2, the equation of state constraints end 
up being nearly independent of the actual value of 7 (see 
Fig. 2 of [13). 

Remarkably, varying the sum of neutrino masses only 
decreases the figure of merit only by 7%. This is mainly 
due to the fact that neutrino free streaming suppresses 
growth of structure in a uniquely scale dependent way, 
with growth on smaller scales (higher k) being more sup- 
pressed than on larger ones, while the effects of dark en- 
ergy, growth index or curvature have no particular scale 
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— none 


1.0 
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---J 
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FIG. 2: Dependence of the dark energy constraints on beyond- 
ACDM parameters: the sum of neutrino masses 'm^, spatial 
curvature 0.^. and the growth index 7. The 68% contours are 
shown for the case where none of these is varied, when each 
parameter is marginalized over individually, and when all of 
them are marginalized over. 



dependence. As a consequence, the correlation coeffi- 
cients between the uncertainties in neutrino mass and in 
dark energy equation of state parameters are small (see, 
e.g.. Table 2 of [H). 

Since this scale dependent suppression of growth makes 
weak lensing power spectra sensitive to the neutrino mass 
sum, we here investigate the constraints, from weak lens- 
ing alone, on the sum of neutrino masses. Figure [3] shows 
the dependence of this constraint, marginalized over all 
other parameters in the WQ-Wa cosmology, on the source 
density of sheared galaxies for fixed survey area (4000 
sq. deg.). Increased source density lowers the uncertainty 
on the shear power spectra, especially on smaller scales 
where the effects of neutrino streaming are more pro- 
nounced. We find changing the sky area scales the con- 
straint by /^Yy . 



C. Redshift Range Contribution 

One of the interesting issues to explore is the red- 
shift range over which weak lensing measurements con- 
tribute significantly to cosmological parameter determi- 
nation. This feeds directly into the survey design trade 
off of a wide area but shallower survey vs a narrower but 
deeper survey. We investigate this in terms of both the 
redshift dependence of the kernel support for the weak 
lensing shear power spectrum (i.e. the lens redshift), and 
in terms of the maximum source galaxy redshift. 

Since the shear power spectrum is an integral over the 
line of sight (Eq. [5]), we can analyze the redshift depen- 
dence of the integrand, or kernel of the projected shear 
power. Recall that z(r]) is the lens redshift. In particular. 
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FIG. 3: Dependence of the constraint on the sum of neutrino 
masses 'm^ on the source density riefi for the 4000 sq. deg. 
survey considered here. 



we expect that for small redshifts the contribution goes 
to zero since -q gets small and for a fixed multipole I the 
scales probed get small and the mass power leading to 
shear decreases. (Alternately, if one considers constant 
fractional distance intervals drj/rj, then the path length rj 
itself gets small.) For high lens redshifts, the lens-source 
distance rjis = r](z, z') gets small and again the integrand 
vanishes. Thus the kernel support peaks at intermedi- 
ate lens redshifts, although the peak does shift to higher 
redshifts as the source distribution, i.e. Zmedj does. (Note 
that the rule of thumb that lensing is strongest at a con- 
formal distance midway between the source and observer 
is more directly relevant to the lensing cross-section than 
to the shear power contribution.) 

Figure d] illustrates the redshift dependence of the 
shear power integrand for several values of multipole t 
and source redshift bin. We show small, medium, and 
large £ and the autopower in bins 1 (0 < z < 0.6), 3 
(1 < z < 1.5), and 5 (2.1 < z < 3.5). Note that the 
peak for (z,) sa 0.3 (1.25, 2.8) is at zio„s < 0.2 (0.4, 
0.7). There is a long tail though to higher lens redshifts. 
Turned around, one could say that to probe gravitational 
lensing in the universe out to z « 0.5 requires sources out 
to z > 2. 

For cosmological constraints, what is of greater impor- 
tance than where the shear power peaks is where lies the 
innate information about the cosmological parameters. 
That is, does information on different cosmologies exhibit 
a similar redshift weighting as the overall shear power in- 
tegrand, or does higher redshift data, for example, show 
greater cosmological leverage? Figure [5] examines this 
issue with respect to the dark energy equation of state 
parameters wq and Wa- 

Considering the finite difference between the kernel 
supports dCi/dz for cosmologies with two different pa- 
rameter values (all other parameters held fixed), e.g. ef- 
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FIG. 4: The contributions of each redshift interval to the 
shear power spectrum, i.e. the kernel support, are plotted vs 
lens redshift to show over what range the power arises. We 
plot t times the shear autopower in source redshift bins 1, 
3, and 5 to show the different multipoles with roughly the 
same amplitude. The interplay of the peak in the mass power 
Pfc and the geometric factor r]is means that high multipoles 
eventually get stronger contributions from smaller redshifts, 
like low multipoles do, while intermediate multipoles probe 
higher redshifts. 

fectively dCe/{dz dwo), the contributions are seen to be 
essentially comparable regarding the peak redshift to the 
shear power itself. For Wa , there is a slight shift in lever- 
age to higher peak redshifts, e.g. the peak for £ = 243, 
bin 3 is z « 0.6 rather than 0.45. The persistence of dark 
energy to higher redshifts through a positive Wa increases 
the cosmology signal at high redshift, while a more nega- 
tive wq diminishes the high redshift dark energy density 
and decreases the signal at high redshift. In both cases, 
however, the rise from small z is pushed to larger z than 
in the power kernel. This makes sense since low redshift 
cosmology is insensitive to the value of wq and even more 
so Wa, and so the leverage is delayed. 

We examine the impact of the survey redshift depth 
on the dark energy equation of state figure of merit in 
Sec. IIVISI where we also investigate the interaction with 
fiducial cosmology. 

D. Source Density and Sky Area 

Another important aspect of survey design concerns 
the number density of source galaxies. The large number 
density of galaxies useful for shear measurements that 
become available with a small point spread function is 
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FIG. 5: The contributions of each redshift interval to the cos- 
mological parameter discrimination ability of the shear power 
spectrum, with respect to wo (top panel) and Wa (bottom 
panel), are plotted vs lens redshift to show from which red- 
shifts the leverage on the parameters arises. The jitter in 
the curves is numerical noise in the derivatives. We multi- 
ply the finite difference between two cosmologies of the shear 
autopower in source redshift bins 1, 3, and 5 hy £ simply to 
plot the different multipoles with roughly the same ampli- 
tude. The redshift dependence of the leverage is similar to, 
or weighted to slightly higher redshifts, than that of the shear 
power itself. 



one of the major advantages of space-based observations 
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[19|. The density rics involves several factors, including 
the PSF characterization and the redshift depth (since 
the number density is projected along the line of sight, 
plus the galaxy size distribution varies with redshift). 

Figure [6] shows that greater source density has a strong 
effect on figure of merit, with almost a linear dependence 
in the range of ricft of interest (the improvement does 
flatten out at higher values, as the shape noise term be- 
comes negligible). This study treats rics as an isolated 
parameter, and does not account for the accompanying 
gain in PSF characterization likely to come from a higher 
density of sources and star calibrations; it also holds the 
survey depth Zmed fixed. The fitting function over the 
range plotted is given by 



FOM(neff) 
FOM(noff = 55) 



O.Oierio 



0.12. 



(27) 
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FIG. 6: The dark energy figure of merit improves significantly 
with source density Ues of measurably sheared galaxies. Halv- 
ing the fiducial source density of 55 measured galaxies per 
square arcminute halves the FOM. 

The total number of lensed sources also scales with the 
sky area, or fraction /sky More sources implies more in- 
formation and so the figure of merit increases with sky 
area. The systematics associated with weak Icnsing, in 
photometric redshifts and in shear measurements, are not 
taken to provide a systematic floor traditionally, i.e. the 
shear power spectrum uncertainty has a factor l//sky 
outside all the error contributions and so continues to 
drop with ever greater area. Whether this is realistic or 
not is an open question. However, priors that are placed 
on marginalized parameters accounting for these system- 
atics do break the linear scaling. 



We examine the dependence of the dark energy param- 
eter estimation and figure of merit on the sky area for the 
particular case of a prior on the photometric redshift bias 
of 10~^. This should be a fairly generic representative 
of the behavior among the systematics. For relatively 
small survey area, the statistical leverage of the survey 
data itself provides the limit on cosmological parameter 
estimation, from fitting for all the systematics, and the 
redshift bias prior is not so important. With large sky 
area, however, the other contributions are well enough 
determined that they become of the same order as the 
prior on Zbias and so the prior plays an important role. 
See, for example. Fig. 7 of The result is that the 

cosmological parameter estimation (FOM) slows down in 
its improvement with /sky, although it does not reach a 
ceiling. 

Figure [7] illustrates the change in scaling. Note that 
even for a survey area of 1000 deg^ the improvement goes 
more slowly than linearly with /sky In the cosmological 
parameter determination there is a noticeable bend in 
the scaling between 4000 and 10000 deg^, after which 
the slope is shallower. A fitting function is 



FOM(area) 



FOM(4000deg ) V4000 



A 



l-0.316(A/4000)" 



(28) 



where sky area A is in units of square degrees. 




10 20 30 

Sky Area (/lOOO sq. deg.) 

FIG. 7: In the presence of priors on systematics (or a system- 
atic fioor), the cosmological parameter estimation does not 
improve linearly with sky area, but somewhat more slowly. 
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E. Effect of Fiducial Model 

One aspect we do not have control over in the survey 
design is the true cosmological model. In the projections 
for how well a survey will constrain the cosmological pa- 
rameters, the true values of the parameters can have an 
impact. Here we examine not only the extent of the influ- 
ence on the parameter estimation but whether one would 
be led to optimize survey design differently depending on 
what one believed the true cosmology was. 

In general, one expects that as the dark energy persists 
to higher redshift, e.g. if its equation of state is closer to 
zero, its effects should be more visible in the shear power 
spectrum and it should be better constrained. Figure [8] 
shows that this is indeed the case. For a shift from w — 
— 1 to less negative values, the area figure of merit in the 
dark energy equation of state WQ-Wa plane increases by 



FOM(w) 
FOM(w; = -1) 



l«6(l + w). 



(29) 



agree with w = —0.9 to better than 4%, and likewise for 
Wq ~ —0.8, Wa = 0.18 and w = —0.86. For fiducial mod- 
els with w more negative than —1, the change in FOM 
has less strong dependence since there is little contribu- 
tion from the higher redshift bins regardless of the exact 
value of w < —1. The scatter in the figure curves near 
If = —1 is a numerical artifact from interaction of the 
small step size (so as not to cross w = —1) with the dark 
energy perturbation evolution. 

These results show that weak lensing surveys may 
have more probative power than expected from taking 
a ACDM fiducial cosmology. However, we have no con- 
trol over what the true cosmology is, so what we do need 
to check is how sensitive conclusions about survey de- 
sign are to the fiducial model. We would not necessarily 
want to choose survey requirements that are near optimal 
for one cosmology but inefficient for another. Figure [9] 
demonstrates that the survey requirements, with regards 
to the survey depth ^mod at least, are robust with respect 
to fiducial cosmology. 



That is, 
1.6. 



a shift Aw = 0.1 boosts the FOM by a factor 




0.8 



FIG. 8: Cosmological parameter estimation improves as the 
effective dark energy equation of state becomes less negative. 
The area figure of merit (FOM) improves by a factor 2 for 
a shift of fiducial cosmology from u; = — 1 to lo = —0.84. 
(Jitter in the curves near u; = — 1 is an artifact of adjusting 
numerical step sizes to prevent crossing of w = —1.) 

Note that here w denotes the effective constant dark 
energy equation of state, defined for the dynamical equa- 
tion of state case using the rule of thumb that the effec- 
tive value w = Wq + Wa/3. We have verified that this 
holds true, e.g. the results for wq = —0.97, Wa = 0.21 




FIG. 9: The dark energy figure of merit improves significantly 
with survey depth, up to Zmed ~ 1, and then nearly saturates. 
This example of survey design optimization is robust against 
assumption of the fiducial cosmology, with the dependence on 
Zmcd being similar for w — —0.9 to —1.1. 

The first important result is that the survey depth Zmcd 
is a crucial factor in the cosmological power of the survey. 
A survey whose median source depth is ^mod = 1 has a 
40% improvement in FOM over one with ^mod = 0.7, all 
else held fixed. However, going deeper than z„icd ~ 1.0 
does not lead to any continued substantial gain. Thus, 
^mcd ~ 1.0 is near optimal for survey design. The second 
important point is that this optimum is not sensitive to 
the (uncontrollable) true cosmology. Models with effec- 
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tive dark energy equations of state between w = —0.9 
and —1.1 all follow this rule, with the FOM as a function 
of Zmod (normalized to the ^mod = 1 case) differing by less 
than 3% between them. An excellent fit for the FOM as 
a function of survey depth over the range plotted is 

FOM(z,„ed = l) \ 0.2 J ^ ' 

V. CONCLUSIONS 

Weak gravitational lensing can be a powerful probe 
into the cosmological model, testing both the expansion 
history and growth of large scale structure. This sen- 
sitivity also requires that a broad array of physics be 
included when calculating the weak lensing shear power 
spectrum and its constraints on cosmological parameters. 
We include the effects of dynamical dark energy, spatial 
curvature, neutrino masses, and gravitational modifica- 
tions into a weak lensing Fisher code, together with sys- 
tematics from photometric redshift measurements, shear 
measurements, and nonlinearities in the mass power spec- 
trum. 

This flexibility allows us to carry out studies analyz- 
ing the effects of each element of an expanded parameter 
space, of systematic uncertainties, and of survey charac- 
teristics on the ability for weak lensing to probe cosmol- 
ogy. We show that allowing growth of structure some 
independence from the expansion history, through the 
gravitational growth index 7, i.e. relaxing the dictum of 
general relativity, enlarges the dark energy equation of 
state WQ-Wa contour area by a factor 2, but could lead 
the detection of new aspects of gravity. We find that in- 
cluding a systematic uncertainty in the nonlinear regime 



behavior of the matter power spectrum has little effect - 
if the k dependence of the modification is known; ignor- 
ing the uncertainty (or getting the dependence wrong) 
can lead to substantial bias in cosmology. 

Systematics control in photometric redshifts and shear 
measurement are, as is well known, critical. The cumu- 
lative effect of these uncertainties is especially powerful. 
The systematics feed strongly into the survey design, for 
example an improvement in shear measurement allowing 
a higher effective number of galaxies yields a sub- 
stantial improvement in the dark energy figure of merit. 

We have presented a number of illustrative survey char- 
acteristic analyses, giving accurate fitting formulas for 
the effect on the dark energy figure of merit with source 
density, survey depth, and sky area. An interesting point 
is that the fiducial cosmology can have a substantial im- 
pact on the figure of merit, but fortunately the behavior 
with survey characteristics (at least as far as depth) is 
nearly independent of fiducial model. Thus survey de- 
sign can focus on systematics, within the control of ex- 
perimental, simulation, and theory efforts. Given these 
efforts, weak lensing promises a cornucopia of knowledge 
about cosmology, gravity, and fundamental physics. 
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